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^ ■ Tip-driven growth processes underhe the development of many complex biological organisms. To date, 

Cn I tip-driven growth processes have been modeled as an elongating path or series of segments without taking 

into account lateral expansion during elongation. Instead, models of growth often introduce an explicit 
thickness by expanding the area around the completed elongated path. Modeling expansion in this way 
rS^ • can lead to contradictions in the physical plausibility of the resulting surface and to uncertainty about 

M-i| how the object reached certain regions of space. Here, we introduce "fiber walks" as a self-avoiding 

random walk model for tip-driven growth processes that includes lateral expansion. In 2D, the fiber walk 
takes place on a square lattice and the space occupied by the fiber is modeled as a lateral contraction 
of the lattice. This contraction influences the possible follow-up steps of the fiber walk. The boundary 

fj I of the area consumed by the contraction is derived as the dual of the lattice faces adjacent to the fiber. 

We show that fiber walks generate a fiber that have well-defined curvatures, enable the identification of 

^L^l the process underlying the occupancy of physical space. Hence, fiber walks provide a base from which to 

(~| ■ model both the extension and expansion of physical biological objects with finite thickness. 



Author Summary 

We consider the growth dynamics (of a part) of an organism whose structure can be represented, as a 
^S| ' single curve. This scenario is inspired by the growth of tap-roots in plants. As a result, we introduce a 

^O , new type of random walk - the fiber walk - that can simulate the growth of thickened objects. We show 

that the random growth of a single extended organism is affected by geometric constraints: it alters where 
^^ , the object can grow, how it grows, and the likelihood that the object is blocked due to prior expansion. 

^^ ' The fiber walk model leads to novel findings in the relationship between growth and form, e.g., the 

fiber walk expansion blocks off regions in space that can never be reached during the growth process, 
but would have been accessible when considering a growth process of "line-like' objects. Compared to 
existing random walk models, the surface associated with the curve is both geometrically and physically 
plausible at all times during the growth process. AUtogether, fiber walks provide a flexible methodology 
}J] , for the physical syudy of organismal growth and form. 

Introduction 

The growth and development of biological organisms involves processes at multiple scales, including regu- 
lation of the timing and location of specific developmental programs. Despite their structural complexity, 
many prototypical processes involved in development have been identified. One such process is the notion 
of tip-driven growth. For example, in the case of plants, root elongation is driven by gradients in the 
growth- hormone auxin that enable elongation at tips ^IfiSj. Above-ground plant and tree structure is 
tip-driven [4], again by gradients in auxin production that enable branch extension at the apical meris- 
tem. Similarly, the long and branching filamentous structure of fungi (hypha) extend from the tips via a 
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mechanism of cell division and physical flow of materials through the hyphae [5] . The fiber walk intro- 
duced in this paper parallels the apical tip-driven growth of individual branches in many organisms due 
to auxin accumulation within an apical meristem [6,- 
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Figure 1. Example of contradictions arising from de-coupling elongation and expansion of a SAW. (A) 
SAW of 16 steps beginning at the origin (yellow triangle) with current position marked via the green 
triangle. (B) Thickened SAW - note the curvature singularity of the resulting surface. (C) SAW 
thickened by 2/3 of an edge length - multiple regions in space are "associated" with different steps, 
leading to an identifiability problem, which holds for all thickening greater than or equal to 1/2 of an 
edge length. 



A number of models have been proposed to examine the growth of extended structures. These vary in 
complexity. For example, perhaps the simplest, from a parameterization perspective, is the self-avoiding 
random walk (SAW) 7_. A SAW is a sequence of steps along the line segments of a grid edges that 
connect the locations of crossing grid-lines vertices. Such a grid is a square lattice in 1? and the SAW 
satisfies the constraint that no vertex is visited twice and all walks are equally likely. 

In SAW and derived models, the elongation and expansion steps are de-coupled. The consequences of 
such de-coupling is that the resulting structure may be non-physical. In Fig[T]we show a SAW thickened 
to 1/3 and 2/3 of an edge length. In both cases, the SAW can be associated with a boundary which is 
exactly the expansion distance away from any point along the edges in the SAW. This derived boundary 
has singular points. Singular points are associated with undefined curvature, which are problematic from 
both geometric and biophysical perspectives. Next, in the case of edge length expansion exceeding 1/2 of 
the edge length, the walk may intersect with itself, thus violating the manifold property that no location 
can be occupied by two surface points. Hence, de-coupling elongation and growth raises the possibility of 
models generating an elongated object that is feasible only in the limit of infinitesimal spatial expansion. 
This violation of the manifold property leads often to difficulties in current reconstruction practice, 8,9 . 

Inspired by the growth of plant roots in particular and by growth processes in biology in general, we 
introduce here the fiber walk: a SAW model that includes a notion of symmetric expansion perpendicular 
to the direction of the path elongation. In an intuitive view a lattice defines a space in which the fiber 
walk grows a fiber along the edges of the lattice. Notably, the interaction of the fiber walk with the 
lattice is two- fold: 1) the elongation of the walk is achieved by stepping to unvisited adjacent vertices on 
the lattice and 2) the spatial expansion is accomplished by contracting certain unvisited vertices adjacent 
to a visited vertex from the side. In comparison to the ordinary SAW, our new contraction allows us 
to define a discrete approximation of the area that is reserved for the fiber to expand. This expansion 
implies a boundary, which is a curve in 2D or a surface in 3D. This defined area gives insight into regions 
of the lattice where no elongation and/or no expansion can occur because of self-avoidance induced by 
the elongation and expansion of a growing object. As we show, the fiber walk resolves the contradictions 
of de-coupled growth models while creating a base for future studies of the relationship between growth 
and form in biology, including branching structures. 



The Fiber Walk 

This section introduces the formal notion of the fiber in n-dimensions, whose process is a walk over the 
edges of a lattice. Our fiber walk includes the notion of a SAW. First we introduce the lattice as a graph 
consisting of vertices connected by edges. Secondly, we define the fiber as a subset of all possible graphs 
on the a lattice whose vertices have at most 2 incident edges. The fiber is characterized by the presence 
of a stopping configuration on the lattice that we define. After introducing the two basic notions of 
lattice and fiber walk, we define fiber walks as the process that creates the fiber as sequence of steps 
along lattice edges. The interaction between the fiber walk and the lattice extends the class of SAW's. 
The extension adds the local contraction of the lattice around the last reached position, which in turn is 
the space consuming expansion of the fiber walk. We illustrate the notions with 2D examples. 
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Figure 2. Example of a growing fiber (green) and an self- avoiding edge (purple) on a lattice as defined 
in Def. [Ij Vertex positions are assigned to each vertex. 



The Lattice 

In 2D a lattice can be thought of infinitely many adjacent squares each of which shares one side and two 
corner points with each adjacent square. These corner points are vertices and the shared sides are edges 
of a graph. In general a lattice, on which the fiber walk grows a fiber, is a hypercube graph L(V, E) 
of infinite size consisting of set of vertices V and a set of edges E. Such a lattice can be derived by 
taking the graph Cartesian product of n one-dimensional lattices at infinite length in Z , where n denotes 
the dimension of the targeted lattice 10 . Beside the notion of dimension, we also need to define an 
embedding for the lattice, which we simply obtain by assigning unit edge length to every edge of the 
lattice. Notably, a position can be assigned to every vertex, such that the vertex positions are denoted 
as n-tuples in Z" (see Figl^]) . The assignment of positions to each vertex will enable us later to define 
contractible edges. 

Definition 1 Let L be a lattice in Z". Fix vq as the origin and assign unit length to all edges, then the 
position of each vertex is the n-tuple associated with Z" . 

The fiber 

We aim at defining the fiber as a non-branching and loop free graph F{V, E) consisting of set of vertices 
V and a set of edges E, [IT], growing on a lattice whose vertices are enumerated in Z". The fiber starts 
at a chosen vertex, the so-called origin. We start with a narrow definition of a path graph directed from 
the origin via incident edges, that excludes the trivial case of the 'empty' path. Furthermore, we restrict 
ourselves in the context of this paper to only one fiber on the lattice. 

Definition 2 Let G be a path graph with m -\- l,m > 1 vertices, such that exactly two vertices have one 
incident edge and m — 1 vertices have exactly 2 incident edges. Fix vq and Vm to be the vertices with one 
incident edge. A path graph is then called directed by the sequence {vo,vi, ...,Vm-i,Vm} 
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Figure 3. Example of a fiber walk on a 2D lattice. Elongation: The first step is chosen randomly 
between the edges incident to the origin vq (compare Fig|4]). Here the chosen edge to reach vi is green. 
Selection of vertices incident to the walk from the side: The vertices selected to be merged with vi are 
shown in grey. Contraction of the selected vertices: The result after merging the selected vertices to vi. 
Another step: The second step on the lattice to reach V2 using the same color schema as before. 

In the following two definitions we introduce the fiber as a subset of possible path graphs. First we 
define self-avoiding edges to limit ourself to cycle-free path graphs and as a measurable characteristic of 
self-avoidance of the walk. Secondly, we use self-avoiding edges to obtain a stopping criteria of the fiber. 

Definition 3 Let L be a lattice and G(V, E) be a directed path graph with m vertices on the lattice. An 
edge e £ L incident to two vertices of G is called a self-avoiding edge if e ^ G. 

Given self-avoiding edges we can define the stopping configuration of a fiber as a vertex that is only 
incident to self-avoiding edges. 

Definition 4 Let L be a lattice. Any directed path graph G — {vo^vi, ...,v„i} with vertices v being a 
subgraph of L is called a stopped fiber, if and only if all incident edges of Vm are self-avoiding edges. 

Note here, that Def|4]only specifies the stopping criterion at v„i and that implicitly the fiber is incident 
to itself, because a fiber contains no vertex with more than 2 incident edges (see Defj^]) . 



The fiber walk 

Until now, we have defined a static view of the fiber, that allows us to identify a fiber on the lattice. 
In the following we describe the process that constructs the fiber. The process of constructing a fiber is 
called a fiber walk and is defined as a sequence of constrained random steps along edges, which can be 
associated with a discrete time step in some applications (e.g. growth of an individual root). The fiber 
walk is a sequence of steps, each from a vertex Vi to an adjacent vertex w^+i via an edge, that is taken by 
an equally random choice of all incident edges to vt to reach ViJ^i. Consequently, at the newly reached 
vertex Vi+i a new random choice is executed. The choice is taken under the constraint that the newly 
chosen edge does not connect to a previously visited vertex. The excluded choices prevent the fiber walk 
from walking back on itself and forming loops. The process is depicted in Fig|4l 

Definition 5 Let Vi be a vertex that belongs to a vertex sequence of 2 or more vertices generated by a 
fiber walk. One step of a fiber walk from a vertex Vi to a vertex w^+i along an edge of a lattice is the equal 
random choice among edges incident to u,; that are not incident to any Vi^k, 1 < k < i. 

Def m -[5] are equivalent to the notion of a growing self-avoiding random walk (FigH]) on a lattice |12) . 
It is sufficient for our purpose to investigate the growing SAW, because its configurations are possible 
configurations of the traditional SAW used to study the excluded volume effect. 
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Figure 4. A step (green) on a lattice and the possible follow-up steps (dotted) 

Fiber walk contraction of the lattice 

The need to explicitly define a contraction as part of the walk is motivated by an identifiability problem 
arising when expansion a walk after elongation is complete (see FiglT|). Here, the expansion of the fiber 
walk causes a contraction on the lattice associated with an elongation step. The contraction of the lattice 
at each step is performed as a set of merges between vertices adjacent to the last reached vertex on the 
lattice. A merge is the union of two adjacent vertices Vi and Vj, where Vi inherits all edges of Vj. A 
complete contraction at Vi is given as a set of merging operations driven by an edge labeling. We first 
define the merging operation of two vertices. As a second step we give the selection criteria to select 
vertices to be merged under the constraint that only edges incident to the walk from the side at the last 
reached vertex are considered for contraction. We denote the merge between two vertices as ©. The 
process is depicted in FiglHl 

Definition 6 Let Vi at position pi be the last vertex reached on the lattice of a fiber walk and let Vj be a 
vertex connected by an non-self-avoiding edge e to Vi. The merge Vi © Vj results in Vi inheriting all edges 
incident to Vj except for e, that is discarded. 

DeflH] describes how the fiber walk interacts with the lattice and implies that the step driving the 
elongation defines the position on the lattice. We recall DefO that defines each step on the lattice as a 
random choice. Notably the choice of the direction can create self-avoiding edges, connecting to previously 
reached vertices. Therefore, we define an initial labeling based on the vertex positions introduced in DeflTJ 
for the initial lattice. After that we give the mechanism to recalculate labels after each contraction step. 
The refined labeling will allow us to limit contraction to edges incident to the walk from the side. 

We assign a bi-directional labeling to the initial lattice edges, implying that each vertex of the lattice 
has the same amount of incoming and outgoing edges. 

Definition 7 Let vi and V2 be two vertices on a lattice with positions p\ and p2 and let both vertices 
be connected by an edge. Two directed edges between vi and V2 exist with labels ±(pi — P2) giving the 
direction of the edge. 

For example, the labels for all SD-directions are: (±1,0,0), (0, ±1, 0),(0, 0,±1). Defl7]is limited to the 
initial labels of the uncontracted lattice. 

The recalculation of the labels considers the merge of two vertices Vi © Vj and is denoted as ffl. Let li 
be the label belonging to the edge ei connecting Vi to Vj. Let the label I2 be an edge 62 7^ ei incident to 
Vj . The label of a newly created edge is then calculated for each entry a in the Cartesian n-tuple of the 
label, such that: 

(1) 
(2) 
(3) 

For example, if li = (1,0, —1) and I2 = (1, 1, —1), then li (B I2 = (1, 1,-1). Recall, a fiber walk is a 
forward-driven process from the a vertex Vi to a vertex w^+i, such that contraction can only occur with 
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edges that are incident to the walk in the current growth direction. Hence, edges entering from the side 
can be identified by an entry in the n-tuple of their label and their length can be directly calculated. 
The selection of edges to contract is therefore given by two rules: A contraction is performed as a set of 
merging operations with vertices having an incident outgoing edge e to u, iff: 

1. e has a outgoing label from w, that differs in its label in more than one value with the last edge added 
to the walk. 

2. e is not an self- avoiding edge 




Figure 5. Boundary of the fiber walk (blue) derived from the face dual of the lattice (red) in 2D for 
the example given in FigO Dotted line segments denote non-unique edges, which do not belong to the 
boundary. 



The expansion of the fiber walk boundary 

In the following we note that a face is bounded by edges that connect vertices. Our central interest is to 
define the spatial expansion of the fiber walk as represented by a boundary. We introduce the fiber walk 
boundary for simplicity in 2D, but the principle extends naturally to higher dimensions. Based on DefEJ 
we can distinguish three kinds of faces sharing an edge or a vertex with the walk: 1.) self-avoiding faces, 
whose bounding edge set contains only edges that are self- avoiding or belonging to the fiber walk, 2.) 
non-self-avoiding faces, having no self-avoidend edges in their bounding edge set or belong to the fiber 
walk and 3.) mixed faces, bounded by edges that are self-avoiding, not self-avoiding or belonging to the 
fiber walk. 

We first introduce the basic definition of the boundary by assuming that no self-avoiding edges are 
present on the lattice. In a second step we will extend this notion to self-avoiding edges. 

FigO shows the derived boundary of the fiber walk on the minimal example given before in FiglS] 

Definition 8 Let L be a lattice, then the dual lattice Ld of L has vertices each of which corresponds to 
a face of L and each of whose faces corresponds to a vertex of L. 

Connsequently connectivity of dual vertices follows directly from the adjacency of the original faces 

Definition 9 Let Vd be the set of vertices v of a dual lattice. Each vertex Vi is connected the vertex w^+i 
if the corresponding dual faces are adjacent. 

As a final component of our setting we extend Def|8]to achieve validity in the presence of self- avoiding 
edges based on an intermediate lattice. The intermediate lattice defines the faces adjacent to the walk 
by placing vertices at a certain fraction of the edge length of incident edges to the walk. This essential 
definition, as we see later on, is depicted in FiglB] 

Definition 10 Let L be lattice and F be a fiber on the lattice. Additionally, let E be the set of lattice 
edges e^ incident to the walk and belonging to the same face and c be the number of contractions causing 
Ci to exist. The intermediate lattice has vertices at maximal distance d{ei) = 2x~c '^^ong e^. 
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Figure 6. (left) the original Lattice with grey edges and orange vertices. The walk is shown in green 
and self-avoidend edges are shown in purple. The black line denotes the half-edge length of edges 
incident to the walk, (right) the intermediate lattice consisting of the black half edge line, the grey 
edges incident to the green walk. Small orange vertices are placed at half-edge distance while the bigger 
orange vertices are the original vertices belonging to the walk. 

Def[TO] defines faces that are incident to consecutive walk steps, such that DefE] and DefE] are appli- 
cable. Furthermore, Def llOl defines the limit case of maximal fiber expansion in which vertices may be 
placed at the same location. The factor c in DeflTUl is needed to compensate for self-avoiding edges of 
length 3, created by two contractions. This special case of unvisited edges on the original lattice corre- 
sponds to regions where an expanding walk can not walk in and expand. We want to stress here, that our 
fiber walk reserves certain regions on the lattice as self-avoidened edges of length 3, but reservation does 
not imply necessarily expansion of the walk. The boundary as defined here contains right angles only in 
the case of an uncontracted lattice, in which all edges are perpendicular to each other. Note, that other 
lattice types may inhibit a right angle on the walk. The given implicit boundary has angles different from 
90° at each boundary vertex where the walk makes a right angle turn, because the faces adjacent to the 
right angle edges do not comply with the self-duality due to the contraction given in DeffTOl We finalize 
this section with actual computations of the introduced fiber walk in 2D and 3D (FiglT]). An example of 
a computed SAW and the fiber walk in 2D and 3D is given in FiglT] FigjT] (a) and (b) show a SAW on 
a lattice including the self-avoiding edges. In comparison, FiglT] (c) and (d) show a fiber walk and the 
lattice resulting from the contraction. Note here, that we generated the lattice only locally around the 
current vertex to avoid memory problems. 

Results 



Fiber boundary does not have curvature singularities 

A fiber walk generates a fiber whose surface is free of curvature singularities, because the contracted 
lattice is not dual to itself. In the case of an SAW, curvature singularities emerge when expansion after 
elongation is complete. For example, in Fig. IB, two adjacent steps at right angles to each other will 
always be associated with a curvature singularity. In contrast, a fiber walk boundary, as seen in FiglS]has 
varying curvature along the entire boundary, even when the walk performs a right angle turn. In all our 
simulations we found two possible lattice configurations at locations where the fiber walk makes a right 
angle turn. The boundary inhibits right angles at the these locations because of the contraction at each 
step given in Dcfl6] We remind the reader that right angles on the boundary can only exist if the lattice 
faces adjacent to the fiber walk are dual to itself. As defined before, self-duality can only occur when the 
fiber makes a right angle turn. Figl5J\ shows locations where the fiber walk forms a right angle. We can 
characterize the two configurations by describing the lattice configuration locally at the side of the right 
angle and assume that no other configuration exists at a right angle turn. Let Vr be a vertex at which 
the two incident walk edges ei and 62 connect to vertices vi and V2 respectively. 




Figure 7. (A) A growing SAW in 2D and (B) a Fiber Walk in 2D. All walk edges are colored in green, 
the lattice is shown in (grey) and the self avoiding edges are colored in purple for all images. (C) A 
growing SAW in 2D and (D) a fiber walk 3D. 

Configuration 1 in Fig|SlA_ has an edge connecting vi and V2- Hence, Def[TU] applies for reconstructing 
the boundary, where the intermediate lattice forms triangular faces at the side of the right angle. Hence 
the face with vertices Vr,vi, V2 is not dual to itself. 

Configuration 2 has an edge e incident to Vr, that does not belong to the fiber and connects to 
a vertex that exclusively belongs to the lattice. Hence, Def|8] and Defl9] apply for reconstructing the 
boundary. The edge e assures that Vr has two triangular faces adjacent, which avoid self-duality. 

By definition of the boundary neither of the two configurations forms self-dual faces adjacent to the 
edges forming the right angle, which makes the fiber walk boundary free of singularities. The resulting 
boundary is a discretized collection of line segments, which we refine with a Catmull-Clark subdivision 
scheme [13] to recover an "organic", i.e. curved, shape. The Catmull-Clark scheme obtains a fine B- 
spline approximation of the boundary, which is valid because no right angles have to be approximated. 
Therefore we can observe that the boundary does not have a curvature singularity and appears as an 
organic shape in Fig HP . 



Fiber walks define a region of occupied space 

A fiber is associated with an explicit region on the lattice defined by the face dual of the lattice. Each 
elongation step of a fiber walk extends the fiber and the boundary. Because the region on the lattice is 
well defined for a given fiber, the area occupied by each step can be identified. Note that not all steps 
are the same size in a fiber walk. Indeed, different step sizes correspond to different degrees of space 
occupancy for each step in a fiber walk. The fiber walk contains 3 edge length classes in 2D of length 
l,'\/2 and 2 (SI Fig. 1), which denotes a step forward on uncontracted lattice edges, a diagonal step 
resulting from the contraction of the lattice and a long step over previously contracted lattice edges. 
Each edge length can be interpreted in terms of the distance the walk must cross through the expanded 
region before elongation can occur again. Note here, that step length over the one-dimensional edges 
varies in our fiber walk model because of the lattice contraction. Each forward step initiates at a point 
inside the actual growing object and has to pass through an expanded region. This initiating point may 
find an analogue in the center of an auxin accumulating region in an individual growing branch. We 
consider the particular classes of step lengths to be dependent on the choice of lattice, but the finding of 
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Figure 8. (A)The illustrated fiber walk (green) is shown on a grey lattice, with two locations marked 
where curvature singularities are avoided. (B) The intermediate lattice (orange). (C) The filled 
boundary (brown) smoothed with a B-Spline in 2D 

varying step length in the fiber walk is generic. See SI Fig 2-3 for 3D and 4D results. 

Fiber configurations generated by fiber walks can be constrained by spatial 
expansion 

Fiber walks generate self-avoiding edges. As explained before, self-avoiding edges correspond to directions 
of growth that are inaccessible as a direct result of the spatial expansion of the fiber. There are five length 
classes of self-avoiding edges for a fiber grown on a 2D square lattice: 1, V2, 2, ^/b and 3 (SI Fig.l). See 
SI Fig. 2-3 for 3D and 4D results. Although these length classes vary with lattice-type, the fact that they 
are heterogeneous and can be linked to spatial expansion is generic. In particular, self-avoiding edges of 
length ^/E and 3 represent regions where expanded regions on the lattice are closer to each other than 
needed for the expansion. Hence, spatial expansion alter the possible configurations of a resulting fiber. 
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Figure 9. Comparison of the growing SAW (blue) and fiber walk (red) stopping times. The figure 
shows the computed walk length at termination of 100.000 single walks. 



Spatial expansion inhibits elongation steps of fiber walks 

The stopping time describes the probability of a walk to generate a stopping configuration (compare 
DefH]). We give here the stopping times for the 2D growing SAW and fiber walk (FigEl) as the cumulative 
density function of the walk length (see SI for computation results). For the computation of the 2D 
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stopping times we considered 100,000 walks of a maximum length of 300 steps. For the growing SAW 
the walks stopped after 50 steps with 50% chance. In contrast the fiber walk reached the 50% mark 
at approximately 29 steps. We found evidence that the two distributions are distinct from each other 
by performing a two-sample Kolmogorov-Smirnov test (p < 0.001,1? = 0.25454). The stopping times 
allow the conclusion that fiber walks have on average less steps until termination then their dimensionally 
reduced equivalent. The computation of exact stopping times in 3D and 4D is beyond scope of this paper. 



Mean Square Distance per step 
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Figure 10. Summarized overview of the MSD scaling of the growing SAW (blue) and the fiber walk 
(green). For both, the fit line in the scaling region is shown in red. 



Expansion lets the fiber walk initially reach further 

From the previous paragraph it follows that we need an indication of how the fiber walk expansion affects 
the growth of the fiber. We compare the end-to-end distance per step to evaluate how far the growing 
SAW and the fiber walk are on average away from the origin. The classic property to characterize this 
behavior is the scaling obtained as the average Euclidean mean-square-displacement (MSD) from the 
origin as a function of the number of steps on the lattice |7j. The scaling of the MSD R{N) in relation 
to the number of steps N taken on the lattice can be denoted as: 
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log{N) 

The scaling exponent d for a walk is determined as the slope of the best fit line where the function R{N) 
reaches an asymptotic slope and is called the scaling region. In contrast the region before reaching an 
asymptotic slope is called the transient region. For example, a standard diffusive process, such as the 
ordinary random walk, corresponds to i? = 0.5. In order to estimate the scaling exponent we computed 
1,000 fiber walks and growing SAW's in dimensions 2,3,4 using simple sampling. Trapped walks were 
restarted by removing the last edge until a walk length of 200 steps was reached to assure that all 1,000 
walks represent configurations with in the scaling region. Note that this restarting technique overrides 
stopping configurations that occur naturally. We selected walks of a certain minimum length to obtain 
the scaling of the MSD over the number of steps on the lattice. The exponents for all scalings given in this 
sections are derived from the least-square fit of a line into the scaling region of the averaged MSD values 
per number of steps. We have chosen a change of less than 0.1 in the asymptotic slope as a threshold 
to determine the scaling region from 30 steps onwards and used bootstrapping to evaluate the error of 
the slope. We obtained for i9. Fig HUK left), a value of 0.55±0.05 and 0.55±0.06 for the 2D fiber walk and 
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the 2D growing SAW respectively. For dimensions 3 and 4 we obtained 0.53±0.04 and 0.53±0.03 for the 
fiber walk and the 0.52±0.04 and 0.53±0.03 for the growing SAW respectively. For computation results 
in 3D see SI Fig. 4 and SI Fig. 5 for the 4D computation. Our main observation is that our fiber walk 
diverges further away from the origin before reaching the scaling region than the growing SAW, which is 
visible in the higher absolute MSD value compared to the SAW. This stronger divergence in the transient 
region suggests that expansion makes fiber walks initially more ballistic because of different step length 
resulting from the contraction. 

A measure for the overall expanded area 

In alignment with the scaling behavior of the self-avoiding edges, the number of merges defines the 
amount of occupied space on the lattice. Nevertheless, the unique characteristic of the fiber walk is the 
contraction, which defines the boundary describing the spatial expansion of the fiber walk. Here we 
give the scaling of the number of contractions for the fiber walk as a measure for the for the size of the 
boundary (See SI Fig. 6-8). In 2D this resulted in an scaling exponent of 1.13, whereas for the 3D and 
4D case a similar exponent of 0.98 and 0.98 was obtained, which tells us that not all edges incident from 
the side result in a merge. Fig. 6-8 in the SI show the actual computation results. 

Discussion 

In this paper we introduced the fiber walk, which is a growing SAW including an expansion perpendicular 
to the growth direction. The expansion of the fiber walk is modeled as a local contraction around the 
last vertex reached on the lattice. We have shown that the fiber walk process constitutes a mechanism by 
which physical space is reserved yet does not imply the expansion of the object into this space. We found 
that the expansion lets a walk diverge initially further away from the origin before entering the scaling 
region and that fiber walk results in walks of less steps until termination than the SAW. Self-avoidance 
as proposed in our model causes encapsulated regions on the lattice inhibiting further exploration by the 
fiber unless a branching mechanism is added. 

In expanding the class of SAWs to include spatial expansion, we compute a non-reversing walk model 
that never moves on edges that lead to a previously visited vertex J14j . Such a model was introduced 
by [15] to model elongating growth, though without the spatial expansion introduced here. One benefit 
of modeling the contraction is to model local object thickness larger than the step size. As far as 
we are aware, the connection between self-avoidance and contraction in SAW-derived models has not 
been established in the study of random walks in biology (see |16j for an overview). In particular, the 
correlated random walk on a lattice has been studied to explain solely the elongation and the direction 
in root development on the example of spruce trees [17) . Recently the importance of spatial expansion 
for modeling plants has been highlighted in [18] and L-Systems might be a suitable formalism to express 
the fiber walk mechanisms underlying underlying the emergence of root shapes [19] . Simulations of the 
root growth have been developed that include rich information on root physiology and biology [20[|21[ . 
yet they do not incorporate the physical and geometric considerations introduced here. 

Future work on the fiber walks will be targeted towards the creation and simulation of physical 
networks including the spatial interactions between several fiber walks on one lattice. We expect the fiber 
model to be influential for the description of below ground growing networks. Many of such networks face 
the problem of being below the maximum resolution of sensing technologies in non-laboratory conditions 
(e.g. roots in real soil). The fiber walk enables us to identify different causes of curvatures and angles, 
which might influence tip-driven growth processes. We expect that such localized descriptors or measures 
may serve as null models in efforts to identify adaptive growth rules in complex organisms. Finally, the 
current definition and implement of the fiber walk best describe the growth of sessile objects. However, 
extensions to the model could also include the possibility of movement of edges within the fiber and 
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further coupling of physical mechanisms of growth influenced by surface properties. The method is 
available within the supporting material of this paper and for download on the webpages of the authors. 

Materials 

The results within this paper were produced with the python program provided together with the paper. 
The software is suitable for developers as a template to integrate our fiber walk into their software. An 
end user should follow the requirements provided in the readme file to execute the program. 
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